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Abstract 

We analytically study the one-dimensional Asymmetric Simple Ex- 
clusion Process (ASEP) with open boundaries under sublattice-parallel 
updating scheme. We investigate the stationary state properties of this 
model conditioned on finding a given particle number in the system. Re- 
cent numerical investigations have shown that the model possesses three 
different phases in this case. Using a matrix product method we calculate 
both exact canonical partition function and also density profiles of the 
particles in each phase. Application of the Yang-Lee theory reveals that 
the model undergoes two second-order phase transitions at critical points. 
These results confirm the correctness of our previous numerical studies. 
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1 Introduction 

The study of steady states of non-equilibrium systems has motivated a lot of 
works over the last decades ^ E]- These models exhibit in general shock 
structures in their steady state. The study of microscopic and macroscopic 
structure of shocks in one-dimensional reaction-diffusion models has also been 
one of physicists' interests recently |H1 [7| • The most well known model 
which reveals shocks in the steady state is Asymmetric Simple Exclusion Process 
(ASEP) with open boundaries The ASEP is a model of diffusing identical 

classical particles with hardcore interactions on a one-dimensional lattice. In 
this model particles are injected from left boundary of a lattice of length L with 
rate a provided that the first site is empty. They hop to right in the bulk of that 
lattice with finite rate and leave it from right boundary with rate f3 provided 
that the last site is already occupied. A substantial amount of exact results have 
been obtained for the ASEP which include the phase diagram, stationary state 
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probability distribution function, correlation functions and correlation lengths. 
It is also known that for special tuning of microscopic rates a travelling shock 
with a step-like density profile evolves in the system. The phase diagram of the 
ASEP with random sequential updating scheme (continuous-time) has three dif- 
ferent phases: a low-density phase, a high-density phase and a maximal-current 
phase. The shocks appear on the coexistence line between the low-density and 
high-density phases, however; for the ASEP with open boundaries the shock 
can be anywhere on the lattice with equal probabilities and the resulting profile 
in the steady state is linear. In order to study the dynamics of the shocks in the 
ASEP, the model been considered on a lattice with periodic boundary condition 
in the presence of a second-class particle (often called an impurity) El E3j • 
In this case the density profile of the particles, as seen from the impurity, is a 
step-like function i.e. a discontinuity between a low-density and a high-density 
region with an interface which is sharp on lattice scale. 

Despite the remarkable work on the microscopic structure of shocks and their 
dynamics in continuous-time update ^JEIE], not much has been known about 
that in discrete-time updates. For the ASEP in sublattice-parallel update the 
temporal evolution of a single shock has been studied |7j . By making use of a 
quantum algebra symmetry of the generator of the process it has been shown 
that on the microscopic level the dynamics of the shock can be described by a 
single particle dynamics in the thermodynamic limit. In the present paper we 
study the microscopic structure of the shocks in the ASEP in sublattice-parallel 
dynamics; however, instead of introducing a second-class particle, we consider 
the case where the total number of particles is kept fixed. 

The ASEP under sublattice-parallel update with open boundaries was first in- 
troduced and studied in In this model particles are injected from the left 
boundary with probability a and extracted from the right boundary with prob- 
ability (3. In the bulk of the lattice they hop to the right with unit probability. 
The exact phase diagram and stationary density profile of the particles in this 
case (without conditioning on the particle number) have been obtained in the 
same reference. Later this model was studied using a matrix product method 
and the same results were obtained ^2). However, the case of fixed particle 
number has not been studied yet. 

Our recent numerical investigations show that in the case of fixed particle num- 
ber the phase diagram of the model highly depends on the total density of the 
particles on the lattice p, a and (3 jlHj . For p <\ the phase diagram contains 
a low-density and a shock phase while for p > ^ it has a high-density phase 
and also a shock phase. In both cases the phase transition points (a c , j3 c ) are 
determined by the total density of particles p. The density profile of the par- 
ticles in the steady state have also been studied for both even and odd sites. 
In low-density (high-density) phase p < ^,(3 > (3 C (p > \.ot> a c ) the density 
profile of the particles for both even and odd sites is a constant on the lattice 
except near the left (right) boundary where it becomes an exponentially increas- 
ing (decreasing) function of site number. In the shock phase p < ~, /3 < /3 C or 
p > | , a < a c the density profile of the particles for both even and odd sites is 
an error function. 

In this paper we continue our studies analytically. Using a matrix product 
formalism, in which the steady state weights are written in terms of the expec- 
tation value of product of non-commuting operators, we calculate the canonical 
partition function of the model exactly. The thermodynamic behavior of this 
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function will also be calculated in each of three phases mentioned above. The 
exact expressions for the density profile of the particles in each phase have also 
been calculated. Recently it has been shown that the classical Yang-Lee the- 
ory of equilibrium phase transitions |14j can be applied to the non-equilibrium 
models to study their out of equilibrium phase transitions (for a review see 
|15j). The simplicity of our model allows us to calculate the line of the canon- 
ical partition function zeros in the thermodynamic limit from which we obtain 
the critical points. In analogy with equilibrium statistical mechanics one can 
introduce the pressure (which is obviously not the physical pressure of the par- 
ticles) and study its behavior as a functions of p _1 for fixed a and (3. It turns 
out that resulting curve looks quite similar to the isotherm of a Bose gas when 
Bose-Einstein condensation takes place. 

Our paper is organized as follows: In section 2 we will define the model and cal- 
culate its steady state probability distribution function using a matrix product 
formalism. The canonical partition function of the model will also be calculated. 
In section 3 the density profile of the particles on the lattice will be obtained. 
In section 4 we will apply the Yang-Lee theory to study the phase transitions 
in the model and also study its similarities with a Bose gas. Finally in section 
5 we will discuss the results. 



2 Canonical partition function 

The ASEP with sublattice-parallel update is defined as follows: Classical par- 
ticles hop deterministically on a one-dimensional lattice of length L which is 
assumed to be an even number. The system is updated in two half time steps. 
In the first half time step the leftmost and the rightmost sites and also the 
particles at even sites are updated. They are injected to the leftmost site with 
probability a provided that the first site is empty. The particles at even sites 
will hop to their rightmost neighbors provided that the target sites are empty. 
If the last site of the lattice is occupied it will be extracted with probability 
0. In the second half-time step only the particles at odd sites will be updated 
according to the bulk dynamics mentioned above. These two steps are shown in 
Fig. 2] From now on we assume that the total density of particles on the lattice 
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Figure 1: The ASEP under sublattice-parallel updating scheme. 
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is kept fixed p — in which M is the total number of particles. In order to 
study the steady state properties of this model analytically we will adopt a so 
called Matrix Product Formalism (MPF). In what follows we will briefly review 
the MPF for the ASEP with sublattice-parallel dynamics first introduced in • 
Let us define the occupation number at site i as Tj S {0, 1} where Tj = 1 if it is 
occupied and Tj = if it is empty. Now according to the MPF, the stationary 
state probability for the system to be in a given configuration {t} = {n, • • • , tl} 
with exactly M particles is given by 

l - 

P({r}) cx S(M - J2 Ti)(W\ flid^OvW) (1) 

i=l i=l 

where we have defined 

&i := nD + (1 - n)E , O t := Ti D + (1 - n )E. (2) 

The Dirac S in guarantees the conservation of the total number of particles. 
For the ASEP with sublattice-parallel dynamics the operators (D, E) and (D, E) 
are finite dimensional square matrices and with the vectors |V) and (W\ satisfy 
the following quadratic algebra 

[E, E] = [D,D] = , ED = [E, D] , DE = 

(W\E(l - a) = (W\E , (W\(aE + D) = {W\D (3) 

(1 - 0)D\V) - t)\V) , (E + 0D)\V) - E\V). 

The operators (D,E) and (D,E) stand for the presence of particles and holes 
at odd and even sites respectively. It has also been shown that Q has a two- 
dimensional representation for a ^ j3 given by |12| 

a _/a(l-/9) 0\ p _( a(3 0\ ( 1 - (3 ' 

(4) 

The normalization factor in which will be called the canonical partition 
function of the model, should be obtained using the fact that the steady state 
probability distribution function P{{t}) is normalized i.e. J2{ T } -P({ T }) = 1- ^ 
can now be written as 

Zl,m = J2 <H M - ]J[02i-i0 2i ]\V). (5) 

{r} i=l »=1 

The canonical partition function of the model Zl.m is a key function in calcu- 
lating the physical quantities such as the density profile of the particles on the 
lattice in the steady (which will be discussed in next section), therefore; we will 
first concentrate on calculating this quantity. As we will see it is more easier 
to calculate so called grand canonical partition function of the model which is 
defined as 

h 

2 

Z l(0 = J2(W\ Y[(r 2l ^D + (1 - T 2t ^)E)(T 2 ^D + (1 - t 2i )E)\V) 

{r} i=l 
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= (W\(CC)i\V)= J2t M Z LtM (6) 

M=0 

where £ is the fugacity associated with the particles and also we have defined 
C := £l) + E and C := ££> + E. The total density of particles on the lattice p 
will then be fixed by the fugacity of them £ through the following relation 

Now one can easily calculate the canonical partition function of the system by 
inverting the series in JHJ using 

where C is a contour which encircles the origin anti-clockwise. 

Let us start with calculating the grand canonical partition function of the model 

using |0J and ©. This can easily be done and one finds 

Z L (0 = (W\(CC)%\V) = 4%) + 4%) (9) 
in which we have defined 



_ ga(a-/3)(l-/3) 

I3{1 - a) - £,a{\ - fif 1 
r(2) u ^ _ (3{a - (3){l - a) 



sno = „, ^ (10) 



where «i = £a 2 (/3 + £(l — /3)) and M2 = /3 2 (^o; + (1 — a)). In the thermodynamic 
limit L — > oo one can easily distinguish three different cases associated with 
three different phases: 

• Case Mi > i*2 

In this case we have — Z^\s). Using J7J) one finds £ = gH~wrj^ 

in this phase. For the fugacity to be positive we should have p> \- Now 
the condition u\ > U2 translates to a > 2(1 — p). These two conditions 
determine the boundaries of this phase which will be called High-Density 
phase hereafter. In order to calculate the canonical partition function of 
the model in high-density phase we apply the standard steepest descent 
method to (JHJ). After some calculations we find the following expression 
for the canonical partition function of the model in this phase: 



_ (a - f3)a L+l (3 L - M {2p - l)*- M +i 
JL ' M ~ (a-2(l-p))(l-/?)*- M (2(l-p)) i -^ 



(12) 



As we will see the density profile of the particles in this phase has an 
exponential behavior near the left boundary and is a constant elsewhere. 

Case Mi < M2 

In this case we have Zl(£) ~ Z^(£). By using Q one finds £ = ^"^"j 
in this phase. As before the fugacity should be positive, therefore; we 
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should have p < ~. The condition u\ < 1*2 also translates to j3 > 2p. 
These two conditions i.e. p < \ and (3 > 2p determine the boundaries of 
the current phase. Using the standard steepest descent method we find 
the following expression for the canonical partition function of the model 
in this phase which will be called Low-Density phase hereafter: 

((3-2p)(2 P ) M +i(i-2p)i- M -i ( ' 

As we will see the density profile of the particles in this phase has an 
exponential behavior near the right boundary and is a constant elsewhere. 

• Case i*i = 1*2 

The remaining phase corresponds to the regions p < ^,(3 < 2p and p > 
\tOl < 2(1 — p). The condition 1*1 = u 2 gives £0 = . It can easily 

be seen that both Z^ 1 and Z^ 1 have a pole at this point. One can also 
verify that the density-fugacity relation Q fails in this phase i.e. the 
fugacity cannot fix the total density of particles on the lattice. It is known 
that the physical meaning for this is the fact that shocks appear in this 
phase. This is the reason that we will call it the Shock phase hereafter. 
In order to calculate the canonical partition function of the model in this 
phase we will again apply the steepest descent method. Let us take the 
contour of the integral a circle with radius R around the origin. The saddle 
point associated with zjp will be £1 = 2 (i-p)[i-p) an< ^ ^ or ^ wm ^ e 
£2 = q^i^p) • I n t nc shock phase we have £1 < £0 < £2- First we take 
the contour of the integral R as R = £2 for both Z^ and Zj* . Then one 
can modify the contour of the integral from R = £2 to R — £1 for Z\/. 
Since Z\} has a pole at £0 there is a contribution from the pole which 
can be calculated using the cauchy residue theorem. The result is that 
the residue Z^ at £2 is the largest contribution to the canonical partition 
function of the model and therefore we have 

Z L . M ~ (a - p)a M p L ~ M (\^)i- M . (14) 

1 - p 

In Fig. |3 we have plotted the phase diagram of the model obtained from above 
discussion. One should note that the dynamical rules in the present model are 
exactly the ones in |llj : however, because of conditioning on the density of 
particles we find a different phase diagram. 



3 Density profile of the particles 

In this section we will define and calculate the exact expressions for the density 
profiles of the particles on the lattice in the steady state. As for the partition 
function it is much easier to calculate the unnormalized average particle number 
in the grand canonical ensemble and then translate the results into those in the 
canonical ensemble using the inversion formula 

,w _ 1 r Jf w£"(Q 



wis. ^.//fW 1 - (15) 
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la 2(l-p) 1 a 



Figure 2: The phase diagrams of the model for p < | (left) and p > | (right). 
The small curves show the behaviors of the density profile of particles in each 
phase obtained from numerical calculations. 



The normalized average particle number at site i should then be obtained from 

(Pi) = {Pi) { l\i/ z lm- (16) 

The average particle number on even sites in the grand canonical ensemble is 
defined as 

(p2i)^(0 = (W\(CC) i - 1 C^D(CC)^- i \V) , l<*<\ (17) 
and for odd sites in the same ensemble 

(P2 i -i)^(O = (w\(ccy- l 0c(cc)^- i \v) , i<i<§ (is) 

in which C and C have the same previous definitions. Using the representation 
of the quadratic algebra J2J) given by l@J and after some algebra one can calculate 
(tT7jl and JTSJ to find 

/ \< u )(c\ £a(a-/3)(l-,3) £ , . 

(0 = "fll-a)-^!-/?)" 1 + (19) 
£a/?(l- a)(a-/3) * 

~r 



(l + (e-l)a)(/3(l-a)-fa(l-/3)) 

£a>(a-/3)(l-a)(l-/3) 
(l + ($-l) a )(i3(l- Q )-(a(l-ffl 



i 2 u 1 



and 



(P2i - 1>L (c) : (20) 

+ _ £ 2 a 2 (l-a)(a-/3)(l-/3) mA ^i 



In what follows we will study the thermodynamic behaviors of (jl9(l and H2()|l in 
each of the above mentioned phases i.e. the high-density, low-density and the 
shock phase. 
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The high-density phase (p > |, a > 2(1 — p)) 

In this phase we have u\ > U2, therefore; in the thermodynamic limit the 
dominant terms in (|19fl are the first and the last term, however; we should 
keep both the first and the second terms in (|2C)|> in the thermodynamic 
limit in this phase. Using fl2")) . lfT5|l . (flf))l and by applying the steepest 
descent method after some algebra one finds 

, \ /o n q(1 - oQ(2p- 1) * 

iP2i - l} = (2P " 1} " 2(1 - ^)(1 - p) + 2(1 -p)) 6 C (22) 

in which i = 1, ■ • • , ^ and that we have defined the correlation length 

C- 1 = |l n 2(l-p)(2(l-/3)(l- P )+a(^2(l-p))) | ^ Ag be geen bQth these 

density profiles are exponentially increasing functions of i near the left 
boundary, however; as we go farther from there they remain constant in 
the bulk of the lattice and also near the right boundary. 

The low-density phase (p < \ , (3 > 2p) 

In this phase we have u\ < 1*2, therefore; in the thermodynamic limit we 
should keep the second and the third terms in (|19fl and only the second 
term in (|2"U)l . Using ljT2l . (|T3|l . i|TT)|) and by applying the steepest descent 
method after some algebra one finds 

(P 2 ,) = 2p + M^e-^ (23) 
4(l- a )(l-/?)p 2 

{p2i - l)= pHi-2 P ) e c (24) 

in which i = 1, • ■ • , -j and we have defined the correlation length = 

In 2al3p ^ s ^i^2p) 1 ' >P I ■ ^ can ^ e seen that the density profiles of the parti- 
cles at both even and odd sites are constant near the left boundary and also 
in the bulk of the lattice while they are exponentially increasing functions 
of i near the right boundary. 

The shock phase (p > 5, a < 2(1 — p) and p < \, (3 < 2p) 
As we mentioned above the density-fugacity relations fails in this phase 
and it can be a sign for the existence of shocks in the density profile of 
the particles on the lattice. The density of particles at even sites on the 
left hand side of the shock (low-density region) is piow-even — and on 
the right hand side of the shock (high-density region) is Phigh-even = 1- 
These values for the density of particles at odd sites are pi ow -odd = 
and Phigh-odd = 1 — a respectively. In what follows we will investigate 
the shock width in more details. In order to calculate the density profile 
of the particles in this phase we adopt the following procedure: for large 
system size the density profile of particles on the lattice can be described 
by a continuous function p{x) in which x = ^ an d < x < 1. From l|19(l 
and H2(J|I one obtains 

(P2 l+ 2) i r ) (o-(P2 l ) i r ) (o^uiur t (25) 

(P2 l+ l)^(0 - (P2r-lA U) (0 OC (26) 
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Now it can easily be verified that for the density of particles at even and 
odd sites we have 

Peven /odd{x, = Peven/ odd {C)e L ' F{xA) (27) 

in which we have defined F(x,£) = (^ — x) lniti + .xlnzi2- Using l|15|) and 
by applying the steepest descent method we find 



d_ 

dx 



Peven/oddi x ) — Peven/ odd^ ^ ' (28) 



where G(x) — (F(x,£ : ) — pln^)|^ = j and £o is the saddle point of the 
integral. It can be shown that G{x) has its maximum value at xq — 
2^2-a-0) ' therefore; one can expand it around xq up to the second order 
to find 

^ P even / oddip^) Peven/odd^- 2 ( ) ( ) (29) 

in which 

r n, ) = -2(2-q-/3) 3 

u W /p( a _ 2(1 - p )) + f3((a - l) 2 + (1 - 2p)) + 2ap(l - a) ' 1 j 

By integrating this expression and applying the above mentioned bound- 
ary conditions we find the following expressions for the density profile of 
the particles at even and odd sites in the shock phase 



Peven(x) = ^ + ^«-/( J 1 1 G" (x ) \ (x ~ X )) (31) 



Poddix) = + V^ er/( v §i G "( x °)i^ - ^ ( 32 ) 

where er/(- • ■) is the error function. From this it follows that the position 
of the microscopic shock position fluctuates around its mean value xq in 
a region of size L 1 / 2 . 



4 Yang-Lee zeros and similarities with a Bose 
gas 

In this section we will first study the applicability of the classical Yang-Lee 
theory to predict the phase transitions and their orders in our model. According 
to this theory the zeros of the canonical partition function as a function of one 
its intensive variables lie on a curve which might interest the real positive axis 
of that variable in the thermodynamic limit in a couple of points [14| . In this 
case the system undergoes a phase transition at these points. The order of 
transition can also be obtained by investigating the angle at which the line of 
the partition function zeros intersects the real positive axis. If the angle is 
then n will be the order of transition. By defining the free energy of the system 
as g = limx,M— >oo t 1h^l,m the line of the canonical partition function zeros 
can be obtained using \IQ 

Re(gi-g 2 ) = (33) 
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in which g\ and gi are free energies of the system on the left and right hand 
side of the transition point. One can make use of (|12|l - 114() , l|33l) and also the 
definition of the free energy to obtain the line of partition function zeros for 
transition between low-density and shock phases and also between high-density 
and shock phase. For the transition between high-density and shock phases 
after some algebra one finds 



(u 2 + v 2 ) 1 



(2(1 -P)) 



{{1-u) 2 + v 2 )%-i (2p-l)P-i 



(34) 



in which we have defined a — u + iv. It can readily be verified that this curve 
intersects the real positive a axis at a c = 2(1 — p) at an angle ^; therefore, 
the phase transition is of second-order. On the other hand for the transition 
between low-density and shock phases we find 



(x 2 +y 2 )% 



(35) 



((1 -x) 2 + y 2 )%-i {l-2p)P-i 
in which we have defined (3 = x + iy. This curve also intersects the real positive 
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Figure 3: The line of the canonical partition function zeros 134(1 (left) and 1)3 5(1 
(right) for p = 0.8 (a c = 0.4) and p = 0.4 ((3 C = 0.8) respectively. 



(3 axis at (3 C — 2p at an angle f ; therefore, as before a second-order phase 
transition takes place. In Fig. |31we have plotted the line of the zeros obtained 
from (PU and lpT5|) . 

Let us now study the similarity between the phase transition in our model and 
that in a Bose gas. In order to see this similarity let us define pressure in terms 
of the grand canonical partition function of the model as P = In Zl(£). Using 
the asymptotic behavior of the grand canonical partition function of the model 
in the thermodynamic limit and 10 one can calculate the pressure as a function 
of total density of particles p. This can easily be done and we find 

[ £M 4 ff P>\,a>2{l-p) 
P = IM^?), P<i/3>2p (36) 

U^^fer 1 ), p>\,a<2{l-p) and p< <2p. 

As can be seen the pressure as a function of density remain a constant in the 
shock phase. It is also a decreasing function of - in both low-density and high- 
density phases. This spectacular behavior reminds us the isotherm of a Bose 
gas when a Bose-Einstein condensation takes place. 
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5 Conclusion 



In this paper we studied the ASEP with sublattice-parallel update and open 
boundaries where the mean particle number is kept fixed. Exact calculations 
using the MPF show that the system undergoes two second-order phase tran- 
sitions. The shock phase exists both for p > | and p < | depending on the 
values of a and [3. For p = \ one has only shock-phase regardless of injection 
and extraction probabilities. In the shock phase the density profile of the par- 
ticles is an error function and its center is located at Iq = 2(2- a- p) 
In order to study the shock dynamics one can consider the ASEP with sublattice- 
parallel update on a ring in the presence of a second class particle. Our pro- 
cedure can also be applied to calcuate the shock profiles in the partially ASEP 
(PASEP) with open boundaries under sublattice-parallel update [P7| . 

6 References 

1. T. M. Ligget, Stochastic Infracting Systems: Contact, Voter and Exclusion 
Processes (Springer, New York, 1999). 

2. P. A. Ferrari, "Shocks in one-dimensional processes with a drift" in "Prob- 
ability and Phase Transition", G. Grimmett (eds.), (Dordrecht: Kluwer, 
1994). 

3. G. M. Schutz, "Exactly solvable models for many-body systems far from 
equilibrium" in "Phase transitions and critical phenomena, vol. 19", C. 
Domb and J. Lebowitz (eds.), (Academic Press, London, 2000). 

4. V. Belitsky and G. M. Schutz, El. J. Prob. 7 Paper No.ll 1 (2002). 

5. B Derrida, L Lebowitz, and E R Speer, J. Stat. Phys. 89 135 (1997). 

6. K. Krebs, F. H. Jafarpour, and G. M. Schutz, New Journal of Physics 5 
145.1 (2003). 

7. C. Pigorsch and G. M. Schutz, J. Phys. A: Math. Gen. 33 7919 (2000). 

8. G. M. Schutz and E. Domany, J. Stat. Phys. 72 277 (1993). 

9. B. Derrida, M.R. Evans, V. Hakim and V. Pasquier, J. Phys. A 26 1493 
(1993). 

10. B. Derrida, S.A. Janowsky, J.L. Lebowitz, E.R. Speer, Europhys. Lett. 
22 651 (1993). 

11. G. M. Schutz, Phys. Rev. E 47 4265 (1993). 

12. H. Hinrichscn, J. Phys. A: Math. Gen. 29 3659 (1996). 

13. F. H. Jafarpour, Phys. Lett. A. 326 14 (2004). 

14. C. N. Yang and T. D. Lee, Phys. Rev. 87, 404 (1952); Phys. Rev. 87, 
410 (1952). 

15. R. A. Blythe and M. R. Evans, Brazilian J. Phys. 33 464 (2003). 



11 



16. S. Grossmann and W. Rosenhauer, Z. Phys. 218 437 (1969); S. Gross- 
mann and V. Lehmann, Z. Phys. 218 449 (1969). 

17. A. Honecker and I. Peschel, J. Stat. Phys. 88 319 (1997). 



12 



